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This  paper  presents  a  novel  system  for  production  of  pure  oxygen  based  on  the  integration  of  a  solid  oxide 
fuel  cell  (SOFC)  and  a  solid  oxide  electrolyzer  (SOEC).  In  the  proposed  arrangement,  the  SOFC  provides 
electricity,  heat  and  H20  in  vapour  phase  to  the  SOEC  which  carries  out  the  inverse  reactions  of  the 
SOFC,  that  is  the  separation  of  H20  into  H2  (used  as  a  fuel  for  the  SOFC)  and  02  (representing  the  yield 
of  the  system).  Simulations  carried  out  in  different  operating  conditions  show  that  when  the  integrated 
SOFC-SOEC  device  runs  at  low  current  densities  (less  than  1000  Am-2),  pure  oxygen  can  be  generated 
with  an  electric  consumption  comparable  to  mid-size  cryogenic  air  separation  units,  and  significantly 
lower  than  small  scale  systems  based  on  the  PSA  technology. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Oxygen,  as  a  pure  component  or  in  high  concentration  mixtures, 
is  used  in  many  field  of  the  industry  as  well  as  in  medical  treat¬ 
ments  [1].  Usually  it  is  obtained  from  air  by  means  of  cryogenic 
and  non-cryogenic  air  separation  technologies  [2,3].  The  cryogenic 
systems  can  fractionate  air  in  a  distillation  process  carried  out  at 
temperature  of  about  -160  °C  by  exploiting  the  different  liquefac¬ 
tion  temperatures  of  02  and  N2.  They  may  deliver  large  oxygen  flow 
rate  (over  3000tonO2  per  day  in  a  single  line)  with  purity  close 
to  100%  and,  among  all  the  air  separation  technologies,  they  are 
the  most  cost-effective  and  the  most  efficient  in  terms  of  energy 
consumption  per  unit  of  product,  when  high  flow  rates  of  oxygen 
are  requested.  Specific  consumption  for  “tonnage”  plants  can  reach 
0.2  kWh  per  kg  of  95%  purity  oxygen  [2,4].  On  the  other  hand,  for 
production  lower  than  about  50tonO2day_1,  alternative  technolo¬ 
gies  become  more  economically  advantageous  [5].  This  is  because 
small  size  cryogenic  plants  have  less  efficient  compressors  and  a 
higher  thermal  input  to  the  cold-box.  As  a  matter  of  fact  in  a  small 
scale  plant  the  specific  consumption  can  be  50%  higher  than  that  of 
a  large  production  unit  [6]. 

Non-cryogenic  systems  are  most  commonly  employed  at  the 
lower  end  of  the  output  rate  scale  and  when  high  purity  oxy¬ 
gen  is  not  needed  (typically  less  than  94%).  They  are  based  on 
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PSA  (pressure  swing  adsorption),  VSA  (vacuum  swing  adsorption) 
technologies,  or  hybrid  solutions  between  the  former  ones,  which 
exploit  the  greater  affinity  for  N2  than  for  02  of  specific  adsorbent 
materials  such  as  the  zeolites.  Such  devices  show  energy  consump¬ 
tion  in  the  range  0.7-1  kWh  kg02_1  for  product  rates  from  800 
to  70  kgo2day-1.  [7,8].  Besides,  other  systems  based  on  chemical 
processes  and  membrane  separation  technologies  are  under  test¬ 
ing,  but  presently  do  not  commercially  compete  with  the  processes 
mentioned  so  far. 

In  principle  oxygen  can  also  be  obtained  through  the  electrolysis 
of  water  (H20);  however  the  high  electrical  consumption  of  this 
process  makes  this  option  not  competitive  with  the  state-of-the-art 
technologies.  A  possible  solution  capable  to  improve  the  efficiency 
of  oxygen  production  via  electrolysis  would  imply  using  H2  (the  co¬ 
product  of  the  water  electrolysis)  as  a  fuel  to  generate  the  electricity 
required  by  the  electrolysis  itself.  This  process  can  be  carried  out  by 
the  integration  of  an  engine  with  a  common  water-electrolyzer,  as 
it  is  reported  in  Fig.  1.  The  operating  principles  are  the  following: 

•  the  water-electrolyzer  allows  the  separation  of  H20  into  the  two 
elements,  hydrogen  and  oxygen.  The  first  is  used  as  a  fuel  for  the 
engine,  while  oxygen  is  stored  in  a  tank. 

•  the  engine  uses  hydrogen  given  by  the  water-electrolyzer  as  a 
fuel  and  ambient  air  as  oxidizer,  producing  part  of  the  electricity 
consumed  by  the  water-electrolyzer. 

However  the  amount  of  H2  generated  through  the  electrolysis 
is  less  than  that  required  by  the  engine  to  produce  the  electricity 
necessary  for  water  dissociation.  This  is  due  to  the  losses  in  the 
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Nomenclature 

A  area(m2) 

BP  SOEC  anode  side  back  pressure  (bar) 

E  specific  work  (J  mol-1 ) 

E°  standard  potential  (V) 

F  Faraday’s  constant:  96487  C  mol-1 

AG  Gibbs  free  energy  change  of  the  reaction  (J  molH2  -1 ) 

h  static  enthalpy  (J) 

h°  total  enthalpy  (J) 

AH  enthalpy  change  of  the  reaction  (J  molH2  -1 ) 

/  current  (A) 

I  work  (J) 

N  number  of  moles 

N  mole  flow  rate  (mol  s-1 ) 

p  pressure  (bar) 

R  universal  gas  constant:  8.314J  mol-1  K-1 

SOEC  solid  oxide  electrolyzer 

SOFC  solid  oxide  fuel  cell 

AS  entropy  change  of  the  reaction  (J  moln2  -1  K_1 ) 

T  temperature  (I<) 

v  velocity  (ms-1) 

V  cell  voltage  (V);  volume  (m3) 

W  electric  power  (W) 

Xj  molar  fraction  of  specie  i  (-) 

Subscripts 
aux  auxiliary 

el  electric 

FIN  final 

is  isentropic 

re  real 

rev  reversible 

spec  specific 

tot  total 


components  and  to  the  fact  that  the  energy  required  for  the  sepa¬ 
ration  of  a  mixture  is  greater  than  zero  even  in  an  ideal,  reversible 
process.  Therefore  in  order  to  allow  a  continuous  operation  of 
the  process  as  shown  in  Fig.  1,  it  is  necessary  to  provide  auxiliary 
energy,  either  as  electricity  to  the  electrolyzer  (shown  in  the  figure 
as  Wel  aux)  or  as  a  fuel  integration  for  the  engine. 

Nevertheless,  the  resulting  energy  consumption  is  still  too  high 
to  compete  with  the  traditional  air  separation  technologies.  This 
is  mainly  due  to  the  relatively  low  efficiency  of  the  electrolyzer 
(70-75%,  calculated  as  the  ratio  between  the  LHV-flow  rate  associ- 
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Fig.  1.  Conceptual  scheme  of  the  process  for  oxygen  production  via  electrolysis. 


Temperature,  °C 


Fig.  2.  Thermal  (TAS),  electrical  (AG)  and  total  (AH)  energy  consumption  for  the 
steam  electrolysis  as  a  function  of  temperature. 


ated  to  the  H2  output  and  the  electric  power  consumption)  and  of 
the  engine  (35-40%). 

As  a  matter  of  fact,  typical  electrolyzers  available  on  the  mar¬ 
ket,  operate  the  dissociation  of  H20  in  liquid  phase  (normally  at 
ambient  conditions):  in  this  case  the  variation  of  the  Gibbs  free 
energy  (AG)  for  the  reaction  H20  ->  H2  + 1/202  -  which  represents 
the  minimum  quantity  of  electrical  energy  required  in  the  reaction 
-  is  237  kj  molH2o-1  •  It  is  worth  to  note  that  the  value  of  AG  lowers 
to  228  kj  molH2o-1  if  the  molecule  of  water  is  in  vapour  phase  (at 
ambient  conditions)  and  decreases  as  temperature  is  raised  (for 
instance  it  becomes  190  kJmolH2o-1  at  750  °C).  Therefore,  from 
a  thermodynamic  point  of  view  it  is  advantageous  to  carry  out  a 
steam-electrolysis  at  high  temperature,  although  this  imply,  on  the 
other  hand,  a  progressive  increase  of  the  heat  required  in  the  reac¬ 
tion.  Fig.  2  explains  this  point:  TAS  and  AG,  respectively  represent 
the  thermal  and  electrical  energies  consumed  in  the  reaction,  while 
AH  corresponds  to  the  total  energy  consumption,  in  the  absence  of 
irreversible  losses.  Although  AH  remains  approximately  constant, 
AG  tends  to  decrease  while  TAS  tends  to  increase  as  the  tempera¬ 
ture  is  raised.  This  means  that  with  the  increase  of  temperature,  it 
is  possible  to  replace  part  of  the  electrical  energy  with  heat,  which 
is,  according  to  the  second  law  of  thermodynamics,  a  lower  grade 
energy  form. 

As  a  consequence,  a  significant  improvement  of  the  system 
shown  in  Fig.  1  can  be  achieved  if  the  cogenerated  heat  produced 
by  the  engine  is  used  together  with  the  generated  electricity  in  a 
steam-electrolyzer  operating  at  high  temperature. 

2.  Integrated  SOFC-SOEC  system 

The  conceptual  scheme  described  in  the  previous  paragraph 
finds  a  very  efficient  application  in  an  integrated  SOFC-SOEC  sys¬ 
tem.  The  SOFC  (solid  oxides  fuel  cell)  [9]  is  a  high  temperature 
(>650  °C)  fuel  cell  that  produces  electric  and  thermal  energy  and 
steam  from  hydrogen  and  oxygen.  The  SOEC  (solid  oxide  elec¬ 
trolyzer)  [10]  is  an  electrolyzer  that  carries  out  inverse  process  of 
a  SOFC.  In  Fig.  3  the  working  principles  of  the  two  components  are 
presented.  In  this  configuration  SOFC  and  SOEC  replace  the  engine 
and  the  H20  splitter  in  the  scheme  reported  in  Fig.  1.  A  detailed 
description  of  the  system  layout  and  the  related  energy  balances 
will  be  reported  in  Section  2.2. 

2. 1  Thermodynamic  principles 

As  it  will  be  better  explained  later,  the  energy  consumption 
of  the  SOFC-SOEC  system  is  due  to  the  difference  between  the 
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Fig.  3.  Working  principles  of  solid  oxides  fuel  cells  (SOFC)  and  solid  oxides  electrolyzer  (SOEC). 


operating  voltage  of  SOFC  and  SOEC.  It  is  therefore  helpful  to 
evaluate  these  quantities  starting  from  the  ideal  conditions.  To  this 
purpose,  a  conceptual  scheme  showing  an  integrated  SOFC-SOEC 
system  is  reported  in  Fig.  4.  In  this  configuration  a  H2-FI2O  mixture 
is  continuously  recirculated  between  the  anode  of  the  SOFC  (where 
H2  is  consumed  and  H20  is  produced)  and  the  cathode  of  the  SOEC 
(where  the  FI20  produced  by  the  SOFC  is  reconverted  into  H2 ).  The 
SOFC  cathode  is  fed  with  ambient  air  where  oxygen  is  considered 
at  a  partial  pressure  of  0.21  atm.  The  result  of  the  process  is  the 
production  of  02  at  SOEC  anode.  By  assuming  that  the  molar  flow 
rate  of  air  fed  at  the  SOFC  cathode  and  of  the  recirculated  FI2-FI2O 
stream  is  significantly  higher  than  the  one  of  the  oxygen  produced, 
it  follows  that  the  composition  of  the  two  former  streams  can  be 
considered  approximately  constant  through  the  cells.  It  is  also 
assumed  that  the  system  operates  at  constant  temperature  T  which 
is  uniform  over  the  SOFC  and  SOEC. 

Under  these  assumptions  the  reversible  potential  can  be  directly 
deduced  from  the  Nernst  equation  for  SOFC  and  SOEC: 
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Fig.  4.  Conceptual  layout  of  an  integrated  SOFC-SOEC  system  for  pure  oxygen  pro¬ 
duction. 


4°vFC  =  £°  +  §  In 


'  ^SOFC 
h2 


,SOFC 


'02 


ySOFC 

Ai_r  n 


(2) 


where  E°  is  the  temperature  dependent  standard  potential,  cal¬ 
culated  at  the  standard  pressure  of  1  atm.  Accordingly,  po2  is  the 
partial  pressure  of  the  oxygen  expressed  in  atm. 

Subtracting  (2)  from  (1)  we  obtain  the  SOEC-SOFC  reversible 
voltage  difference: 


i/SOEC  _  i/SOFC 
v  rev  v  rev 


pS°ec\ 

PSo°FCJ 


(3) 


If  the  current  /  circulates  in  the  cells,  the  resulting  power  consump¬ 
tion  of  the  ideal  SOFC-SOEC  system  will  be: 


Wrev  =  I  (Vf>°vEC  -  CvFC)  =  I§  In  W 

Since  the  current  I  can  be  expressed  as  a  function  of  the  mole  flow 
rate  of  the  02  output  as: 


I  =  No24F  (5) 

and  by  combining  Eqs.  (4)  and  (5)  we  obtain  the  specific  energy 
required  for  oxygen  production: 


Frev  — 


wrev  ijV^-V^)  RT  |n 
No,  N0,  4 F  n 


'  nSOEC N 
p02 

rl0K 


(6) 


which  is  equivalent  to  the  ideal  work  required  to  separate  and  make 
available  at  pressure  Pq°ec  the  oxygen  contained  in  the  ambient  air 
at  the  partial  pressure  Po°FC  (assumed  as  0.21  bar)  in  a  isothermal 
process  carried  out  at  temperature  T  (see  Appendix  A). 

With  reference  to  the  value  of  p^EC,  three  cases  are  worth  to  be 
considered: 


a)  PS£EC>PS0°JC 

This  is  the  case  where  oxygen  is  produced  by  the  SOEC  for 
instance  at  ambient  pressure.  This  solution  brings  about  the 
advantage  of  avoiding  any  further  O2  compression;  on  the  other 
hand  the  power  absorbed  by  the  SOFC-SOEC  is  grater  than  zero 
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Fig.  5.  Examples  of  voltage-current  curves  for  SOFC  and  SOEC  operating  at  different 
pressures. 


and  proportional  to  the  operating  temperature  according  to  Eq. 

(6). 

(b)  ps0°2EC=ps0°2¥C 

The  oxygen  produced  by  the  SOEC  has  the  same  pressure  of 
the  oxygen  in  the  air  and  ideally  no  power  is  supplied  to  the 
SOFC-SOEC  system.  Power  is  required  subsequently  to  bring 
the  oxygen  to  ambient  pressure.  This  solution  is  more  efficient 
than  case  (a)  because  the  compression  process  can  be  carried 
out  after  cooling  02  at  ambient  temperature  thus  reducing  the 
power  required  according  to  Eq.  (6). 

fr\  nSOEC  ^  nSOFC 

Po2  ^02 

The  SOFC-SOEC  system  generates  power  which  in  an  ideal 
case  is  equal  to  the  power  that  can  be  obtained  in  an  ideal 
isothermal  expansion  of  02  from  p^ofc  (0.21  bar  in  this  case) 

to  pgfc.  This  is  possible  only  if  an  equal  amount  of  heat  is  pro¬ 
vided  to  the  SOEC  in  order  to  equilibrate  the  energy  balance  of 
the  SOFC-SOEC  system. 

Aiming  to  better  focus  these  aspects,  Fig.  5  shows  in  the  V-I  plane 
a  typical  curve  of  a  SOFC  (operating  at  ambient  pressure)  together 
with  three  curves  of  SOECs  characterized  by  different  anode  side 
back  pressures.  Curve  SOEC  BP  =  1.01  refers  to  the  case  where  oxy¬ 
gen  is  produced  at  ambient  pressure.  It  can  be  seen  that  at  open 


circuit  (J  =  0)  Vfe°EC  is  higher  than  Vrse°EC;  in  this  condition  the  dif¬ 
ference  of  the  two  voltages  can  be  calculated  according  to  Eq.  (3) 
and  yields  about  0.034  V  (at  T=750K).  Curve  BP  =  0.21  and  BP  =  0.1 
refer  to  the  cases  where  back  pressure  is  p  =  0.21  bar  and  p  =  0.1  bar 
respectively.  For  1  =  0  curves  SOFC  and  SOEC  BP  =  0.21  have  the 
same  reversible  potential,  while  a  negative  difference  Vrse°EC  -  Vrse°EC 
results  for  the  case  BP  =  0.1.  On  real  operating  conditions,  when 
the  device  runs  at  />  0,  activation,  ohmic,  and  concentration  losses 
take  place,  causing  a  progressive  increase  of  the  SOEC  voltage  (and 
correspondingly  a  decrease  in  the  SOFC  voltage).  This  leads  to  an 
increase  of  AV  between  the  two  components  thus  lowering  the 
oxygen  production  efficiency. 

2.2.  System  layout 

In  Fig.  6  the  flowsheet  of  the  integrated  SOFC-SOEC  device  is 
shown.  In  this  scheme,  the  SOEC  anode  side  is  fed  (stream  8)  with 
a  H20  rich  stream  coming  from  the  SOFC  anode.  Inside  the  SOEC, 
oxygen  is  separated  and  sent  to  storage  (stream  14).  A  multi-stage 
intercooled  compressor  is  employed  to  bring  the  oxygen  stream  to 
ambient  pressure  when  the  SOEC  anode  operates  at  pressure  lower 
than  atmospheric.  The  stream  at  the  SOEC  anode  side  outlet  (rich  in 
H2,  point  7)  is  addressed  to  the  SOFC  anode  for  power  generation. 
It  should  be  noted  that,  due  to  the  comparatively  high  operating 
temperature  (about  750  °C),  instead  of  a  traditional  blower  an  ejec¬ 
tor  ensures  the  recirculation  of  the  H2-H20  mixture  from  the  SOEC 
cathode  outlet  to  the  SOFC  anode  inlet.  To  this  purpose,  stream  9, 
a  fraction  of  the  stream  at  SOEC  cathode  exit,  is  sent  to  the  heat 
exchanger  where  it  is  cooled  down  to  a  temperature  (stream  10) 
acceptable  for  a  compression  carried  out  by  a  conventional  blower. 
The  heat  exchanger  also  increases  the  system  efficiency  by  preheat¬ 
ing  the  air  flow  at  the  SOFC  inlet  (stream  3)  with  heat  recovered 
from  the  streams  of  depleted  air  (stream  4)  and  oxygen  (stream  12) 
exiting  the  SOFC  and  SOEC,  respectively.  The  figure  also  shows  the 
electricity  and  heat  flow  transferred  from  the  SOFC  to  the  SOEC. 
With  reference  to  the  latest  aspect,  the  “sandwich”  configuration 
proposed  in  Fig.  7  (obtained  by  overlapping  layers  of  SOFC  and  SOEC 
realized  according  to  a  planar  technology  [  11,12  ] )  could  be  a  solution 
to  ensure  an  effective  heat  transfer  between  the  components. 

Moreover  it  is  worth  to  note  that  for  assigned  values  of  cell  active 
area  and  current  density  -  being  therefore  fixed  the  02  yield  -  the 
solution  with  a  high  number  of  SOFC  and  SOEC  cells  of  reduced 
area  could  be  beneficial.  In  fact  this  determines:  (i)  a  reduction  of 
the  current  flowing  across  the  stack  which  in  turn  brings  about  a 


Fig.  6.  Flowsheet  of  an  integrated  SOFC-SOEC  system  for  pure  oxygen  production. 
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Fig.  7.  Possible  arrangement  for  assembling  the  SOFC  and  SOEC  stack. 


reduction  of  the  ohmic  losses  -  proportional  to  the  square  of  the 
current1  -  in  the  circuit  connecting  the  auxiliary  electric  generator 
with  the  SOFC-SOEC  stack  (Fig.  7);  (ii)  lower  heat  losses  due  to  the 
lower  value  of  the  surface  to  volume  ratio  of  the  stack.  In  fact  the 
thickness  of  a  planar  cell  is  typically  much  lower  than  the  values 
of  its  length  and  width:  therefore  by  assembling  the  stack  with  a 
higher  number  of  small  size  cells  it  is  possible  to  obtain  a  signif¬ 
icant  reduction  of  the  area  of  the  bottom  and  upper  surface,  with 
a  comparatively  small  increase  of  the  lateral  stack  surface.  These 
considerations  are  reasonable  within  a  certain  number  of  cells  in 
the  stack  as  the  benefit  overcomes  the  disadvantages  related  to  the 
increasing  costs  for  interconnections  and  fluids  distribution  and  to 
the  reduced  reliability. 

As  far  as  thermal  balance  is  concerned,  Fig.  7  shows  that,  owing 
to  the  endothermic  and  exothermic  nature  of  the  reactions  within 
SOFC  and  SOEC,  respectively,  heat  is  transferred  by  conduction  from 
the  a  SOFC  cell  to  the  two  adjoining  SOECs  in  the  stack,  while  with 
regard  to  the  electrical  connections  it  can  be  noted  that: 

-  SOFCs  and  SOECs  are  connected  in  series,  meaning  that,  in  steady 
state  conditions,  the  current  remains  the  same  through  the  stack. 
(By  definition  the  current  flows  in  the  opposite  direction  with 
respect  to  the  electrons  flow  “ 2e~ "  shown  in  Fig.  7.) 

-  Direct  contact  between  the  SOEC  cathode  and  SOFC  anode  on  one 
side,  and  between  the  SOFC  cathode  and  the  SOEC  anode  on  the 
other  side  simplifies  the  electrical  connection. 

-  Since  in  typical  operating  conditions  the  voltage  delivered  by  the 
SOFC  is  lower  than  that  required  by  the  SOEC  it  is  necessary  to 
supply  the  required  additional  voltage,  by  means  of  an  auxiliary 
generator;  as  already  mentioned  this  component  is  responsi¬ 
ble  of  the  energy  consumption  of  the  device.  Otherwise,  the 
lack  of  potential  can  be  compensated  providing  additional  SOFC 
cell  to  the  stack,  which  would  be  fed  by  a  supplementary  fuel 
source. 

According  to  the  previous  indications,  the  proposed  system 
shows  the  following  advantages: 

•  The  SOEC  can  dissociate  the  water  molecules  at  high  tempera¬ 
ture  and  in  vapour  phase.  Therefore  it  substantially  reduces  the 


1  For  a  fixed  value  of  the  electric  power,  this  implies  a  higher  voltage  provided  by 
the  auxiliary  generator. 


electric  absorption  compared  to  the  ordinary  electrolyzers  that 
operate  at  low  temperature  on  liquid  water. 

•  The  SOFC  generates  electricity  at  very  high  conversion  efficiency 
(50-70%).  The  remaining  fuel  input  is  released  as  high  tempera¬ 
ture  heat  available  for  SOEC  (the  vapour  H20  dissociation  reaction 
is  endothermic,  the  heat  transferred  from  SOFC  to  SOEC  allows  to 
keep  the  system  at  the  correct  operating  temperatures). 

•  The  system  does  not  consume  any  water. 

3.  Simulation  model  and  results 

To  estimate  the  performance  of  the  integrated  SOFC-SOEC  sys¬ 
tem  depicted  in  Fig.  6,  the  main  operating  conditions  and  design 
parameters  have  been  set  and  the  complete  energy  balance  has 
been  calculated  by  means  of  a  computer  simulation.  A  75%  utiliza¬ 
tion  factor2  (for  H2  in  the  SOFC  and  H20  in  the  SOEC  respectively) 
has  been  assumed  by  setting  the  flow  rate  in  the  loop  according  to 
the  generated  02  flow  rate. 

The  electric  performance  of  the  cells  has  been  defined  by  means 
of  the  specific  polarization  curves  shown  in  Fig.  5,  where  the  volt¬ 
age  at  open  circuit  (/  =  0)  is  determined  for  SOFC  and  SOEC  by  Nernst 
equations  (Eqs.  (1)  and  (2))  assuming  a  mean  value  of  0.5  for  the 
molar  fraction  of  H2  and  H20.  The  complete  polarization  curves 
of  SOEC  are  then  obtained  considering  the  cell  overpotential  as 
a  function  of  the  current  density,  taken  from  a  simulation  model 
developed  for  an  IT-SOEC  by  Udagawa  et  al.  [10].  In  case  of  SOFC 
it  is  assumed  that  the  basic  components,  materials  and  electrical 
properties  and  operating  temperature  are  the  same  as  SOEC  and 
accordingly  it  is  reasonable  to  consider  that: 

9yS0Fc  dysoEc 
di 

meaning  that,  for  a  given  current  density  /,  the  electrical  losses 
which  determines  the  change  from  the  Nernst  potential  are  the 
same  for  both  SOFC  and  SOEC. 

The  pressure  at  which  oxygen  is  separated  within  the  SOEC 
is  matter  of  optimization  both  from  thermodynamic  and  tech¬ 
nologic  perspective.  For  an  assigned  final  pressure  of  02  stream 
(14  in  Fig.  7),  in  order  to  minimize  the  specific  energy  consump¬ 
tion  it  is  necessary  to  minimize  the  sum  of  the  power  supplied 
by  the  auxiliary  generator  and  the  one  required  by  the  final  com¬ 
pression  which  respectively  increases  and  decreases  with  02  back 
pressure  on  the  SOEC  anode.  Therefore,  for  a  fixed  value  of  cur¬ 
rent  density  in  the  stack,  the  optimized  oxygen  pressure  depends 
on  SOFC  and  SOEC3  polarization  curves  and  on  the  features  of 
the  final  compression  (namely  number  of  intercooled  stages  and 
efficiency).  According  to  the  calculation  assumptions  reported  in 
Table  1  the  analysis  shows  that  any  pressure  reduction  has  a  bene¬ 
ficial  effect  on  the  system  efficiency.  A  back  pressure  value  of  0.1  bar 
has  been  selected  as  preliminary  assumption  to  contain  the  num¬ 
ber  of  compression  stages  and  the  volume  flow  rate  on  the  SOEC 
anode  side  even  if  a  complete  design  of  the  SOFC-SOEC  device 
should  include  a  mechanical  stress  analysis  aiming  to  assess  the 
electrolyte  resistance  to  the  pressure  difference  applied  between 
anode  and  cathode.  Under  these  assumptions,  the  resulting  values 
of  SOFC  and  SOEC  Nernst  potentials  are  0.955  and  0.939  V,  respec¬ 
tively. 

Finally,  the  ejector’s  driving  flow  pressure  and  mass  flow 
(stream  11)  requested  to  overcome  the  pressure  losses  of  the 


2  The  utilization  factor  is  the  ratio  between  the  reacted  moles  of  a  species  and  the 
number  of  moles  of  that  species  at  the  component  inlet. 

3  In  principle,  stream  back  pressure  on  the  anode  side  may  affect  the  SOEC  polar¬ 
ization  curve.  This  second  order  effect  has  not  been  considered  in  the  present 
discussion. 


P.  Iora,  P  Chiesa  /  Journal  of  Power  Sources  190  (2009)  408-416 


413 


Table  1 

Simulation  assumptions. 


SOFC-SOEC  current  density  (A  rrr2 )  500 

SOFC  operating  voltage  (V)  0.9179 

SOEC  operating  voltage  (V)  0.9760 

Pressure  losses  in  air  piping  (streams  2-5  in  Fig.  6)  (Pa)  2000 

Pressure  losses  in  H2-H20  loop  (streams  6-8  in  Fig.  6)  (Pa)  400 

Minimum  AT  within  the  recuperator  (°C)  50 

Cold  blower: 

Isentropic  efficiency  0.7 

Organic/electric  efficiency  0.8 

Stream  fraction  to  ejector  driving  flow  0.05 

Ejector  circuit  compressor: 

Isentropic  efficiency  0.7 

Organic/electric  efficiency  0.8 

02  intercooled  compressor: 

Number  of  stages  3 

Intercooling  exit  temperature  (°C)  35 

Cooler  pressure  drop  (Pa)  2000 

Isentropic  efficiency  0.75 

Organic/electric  efficiency  0.9 


H2-H20  loop  (streams  6-8  in  Fig.  4)  are  calculated  according  to 
the  analysis  carried  out  in  Appendix  A. 

Based  on  these  assumptions,  simulations  of  the  system  shown 
in  Fig.  6  have  been  carried  out  with  the  ASPEN  PLUS  software. 
Properties  of  each  streams  in  the  relevant  points  of  the  flowsheet 
in  Fig.  6  are  reported  in  Table  2  while  Table  3  reports  the  ther¬ 
mal  balance  summary.  According  to  the  above  defined  assumption 
set,  for  a  current  density  of  500  Am-2  the  energy  consumption  is 
about  0.30 kWh  per  kg  of  generated  02.  It  is  worth  to  note  that 
when  the  SOEC  anode  operates  at  atmospheric  pressure,  the  energy 
consumption  rises  to  0.39  kWhkg02_1,  a  value  comparable  to  the 
performance  achievable  by  mid-size  cryogenic  air  separation  units, 
and  significantly  lower  than  the  small  scale  units  based  on  the  PSA 
technology.  Even  higher  efficiency  can  be  achieved  by  adopting 
high  temperature  SOFC  and  SOEC  that  operate  at  about  1000  °C  (the 
estimations  here  reported  refer  to  intermediate  temperature  cells 
operating  in  the  range  650-800  °C),  or  by  assuming  lower  SOEC 
anode  side  back  pressure. 

On  the  other  hand  an  increase  of  the  current  density  brings  about 
a  higher  02  flow  rate  but  also  a  strong  increase  of  the  electric  con¬ 
sumption  due  to  the  higher  polarization  losses  and  the  reduced 
recuperator  effectiveness.  This  can  be  seen  from  Fig.  8  where  the 
difference  between  the  SOEC  and  SOFC  voltage  profile  (Fig.  8A) 
-  which  is  the  relevant  parameter  for  the  power  consumption  of 
the  SOFC-SOEC  system  -  and  the  specific  consumption  (Fig.  8B) 
are  plotted  as  a  function  of  the  current  density  in  case  of  02  back 
pressure  on  the  SOEC  anode  of  1.01  bar  and  0.1  bar,  respectively.  It 


Table  3 

Main  results  of  the  energy  balance  of  the  integrated  SOFC-SOEC  system. 


Specific  electric  energy  produced  by  the  SOFC  3.075 

(kWhkgo2-1) 

Specific  electric  energy  consumed  by  the  SOEC  3.270 

(kWhkg02_1 ) 

Specific  electricity  consumed  by  air  blower  (kW h  kg02  -1 )  0.00928 

Specific  electricity  consumed  by  the  ejector  circuit  0.00860 

compressor  (kW  h  kg02  _1 ) 

Specific  electricity  consumed  by  the  02  intercooled  0.0885 

compressor  (kW  h  kg02  _1 ) 

Specific  consumption  (kW  h  kg-1  of  pure  02 )  0.3010 

Specific  consumption  (kW  h  Nm-3  of  pure  02 )  0.4297 

Heat  exchanged  between  SOFC  and  SOEC  (kW  h  kg02  )  1 .036 

Heat  exchanged  within  the  recuperator  (kW h  kg02  -1 )  2.329 


Current  density,  A.m'2 


Fig.  8.  Difference  between  the  SOEC  and  SOFC  voltage  profile  (A)  and  the  specific 
consumption  (B)  plotted  as  a  function  of  the  current  density  in  case  of  SOEC  anode 
back  pressure  (BP)  0.1  bar  and  1.01  bar. 


Table  2 

Properties  of  significant  streams  in  the  flowsheet  shown  in  Fig.  6. 


Point 

°C 

bar 

kg  kgo2 1 

kmolkgo2“1 

Ar 

h2 

h2o 

n2 

02 

1 

15.0 

1.013 

11.50 

0.3986 

0.0100 

0.0000 

0.0100 

0.7730 

0.2070 

2 

17.3 

1.033 

11.50 

0.3986 

0.0100 

0.0000 

0.0100 

0.7730 

0.2070 

3 

700.0 

1.028 

11.50 

0.3986 

0.0100 

0.0000 

0.0100 

0.7730 

0.2070 

4 

750.0 

1.018 

10.50 

0.3674 

0.0109 

0.0000 

0.0109 

0.8388 

0.1395 

5 

79.7 

1.013 

10.50 

0.3674 

0.0109 

0.0000 

0.0109 

0.8388 

0.1395 

6 

750.0 

1.013 

0.5162 

0.09896 

0.0000 

0.8000 

0.2000 

0.0000 

0.0000 

7 

737.9 

1.017 

0.5433 

0.1042 

0.0000 

0.8000 

0.2000 

0.0000 

0.0000 

8 

750.0 

1.015 

1.543 

0.1042 

0.0000 

0.2000 

0.8000 

0.0000 

0.0000 

9 

750.0 

1.013 

0.02717 

0.005209 

0.0000 

0.8000 

0.2000 

0.0000 

0.0000 

10 

350.0 

1.008 

0.02717 

0.005209 

0.0000 

0.8000 

0.2000 

0.0000 

0.0000 

11 

503.2 

1.823 

0.02717 

0.005209 

0.0000 

0.8000 

0.2000 

0.0000 

0.0000 

12 

750.0 

0.100 

1.000 

0.03125 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

13 

79.7 

0.098 

1.000 

0.03125 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

14 

35.0 

1.013 

1.000 

0.03125 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

Mass  and  molar  flow  rates  are  referred  to  1  kg  of  02  output.  Compositions  are  given  as  mole  percent  of  the  total  flow. 
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Fig.  A.l.  Conceptual  schematics  of  an  ideal,  isothermal  process  to  make  available  oxygen  at  pressure  pFiN.  from  ambient  air. 


can  be  seen  that,  assuming  for  instance  a  current  density  twice  the 
base  value  (1000  vs.  500  A  m-2 )  and  keeping  the  temperature  of  the 
cells  at  750  °C  by  controlling  the  air  flow  rate  on  the  SOFC  cathode 
side,  the  oxygen  production  doubles,  but  the  specific  energy  absorp¬ 
tion  substantially  increases  from  0.30  to  0.55  kWh  kg-1  of  pure  02 
for  BP  =  0.1  bar,  and  from  0.39  to  0.64  kW  h  kg02  -1  for  BP  =  1.01  bar. 
However  these  results  can  still  be  acceptable  if  compared  with  low 
size  oxygen  production  devices  based  on  PSA  technology. 

As  a  final  remark  it  is  worth  to  emphasize  that,  being  a  typical 
modular  system,  this  device  allows  to  reach  the  desired  02  prod¬ 
uct  output  by  simply  increasing  the  number  of  parallel  modules 
without  any  performance  decrease.  Given  that  no  significant  scale 
effects  can  be  predicted,  small  size  appliances  are  expected  to  be  the 
first  commercial  application  of  this  technology  (output  in  the  range 
ofl-101min-1  of02  typical  of  home  oxygen  therapy  treatments,  or 
alternatively  for  02  production  in  hospitals).  As  far  as  the  size  of  this 
appliance  is  concerned,  a  gross  estimate  of  the  bulkiest  part  of  the 
system  (the  cell  stack)  can  be  done  in  case  of  small  prototype  with  a 
2.5 1  of  02  per  minute  yield.  With  a  1000  A  m-2  current  density  and 
assuming  a  stack  of  20  planar  cells  (10  operating  as  SOFC  and  10 
operating  as  SOEC),  it  results  a  block  with  a  720  cm2  cross  surface 
(approximately  a  27  cm  side  square)  and  about  50  cm  high,  while 
a  recuperator  heat  transfer  surface  of  0.063  m2  can  be  estimated.4 
This  would  result  fully  comparable  to  commercial  systems  for  small 
02  production  used  for  home  oxygen  therapy  treatments. 

4.  Conclusions 

This  paper  described  an  integrated  SOFC-SOEC  device  able  to 
generate  pure  oxygen  with  an  electric  consumption  comparable  to 
mid-size  cryogenic  air  separation  units,  and  significantly  lower  than 
small  scale  systems  based  on  the  PSA  technology. 

This  conclusion  is  the  result  of  an  investigation  that  accounted 
for  the  real  operating  conditions  of  the  components.  In  particular, 
the  performance  of  SOFC  and  SOEC  at  the  design  point  has  been 
assumed  according  to  definite  polarization  curves  of  these  compo¬ 
nents.  The  energy  balance  of  the  complete  system,  shows  an  electric 
specific  consumption  of  about  0.30  kWh  kg-1  of  pure  02  output 
when  the  system  is  operated  at  a  current  density  of  500 Am-2. 
Besides  the  clear  advantages  in  term  of  energy  and  money  sav¬ 
ings,  other  favorable  aspects  can  be  envisaged  over  the  competing 
technologies: 


4  Calculated  by  assuming  an  overall  heat  transfer  coefficient  of  250Wm2K_1, 
average  value  for  a  compact  heat  exchanger  operating  at  ambient  pressure. 


(1)  The  SOFC-SOEC  technology  virtually  produces  100%  pure  oxy¬ 
gen.  Such  a  result  is  not  achievable  in  small-scale  commercial 
system  (typically  based  on  PSA  processes  that  allows  to  attain 
02  purity  in  the  range  90-94%).  For  the  large-scale  system,  the 
multiple  columns  configuration  allows  to  achieve  a  purity  close 
to  100%  but  it  entails  an  increase  of  the  energy  consumption. 

(2)  Based  on  the  fuel  cell  technology,  the  device  is  typical  modular 
system.  It  is  then  possible  to  reach  the  desired  02  product  output 
by  simply  increasing  the  number  of  parallel  modules  with¬ 
out  any  performance  decrease.  Given  that  no  significant  scale 
effects  can  be  predicted,  small  size  appliances  are  expected  to 
be  the  first  commercial  application  of  this  technology  (output 
in  the  range  of  1  -10 1  min-1  of  02  typical  of  home  oxygen  ther¬ 
apy  treatments,  or  alternatively  for  02  production  in  hospitals), 
also  considering  that,  in  the  case,  a  considerable  energy  sav¬ 
ing  would  be  achieved.  From  a  preliminary  design,  the  size  of 
the  device  resulted  fully  comparable  to  commercial  systems  for 
small  02  production  used  for  home  oxygen  therapy  treatments. 

(3)  It  is  also  possible  realizing  a  transportable  appliance  where 
the  energy  requirement  is  derived  from  liquid  fuels  (methanol 
or  even  gasoline)  or  provided  by  hydrogen  stored  in  a  bottle. 
With  reference  to  the  latest  option  and  according  to  the  results 
obtained,  the  device  would  be  able  to  produce  1  Nm3  of  pure 
oxygen  using  0.222  Nm3  of  H2.  Given  that  a  bottle  of  pressur¬ 
ized  gas  can  hold  approximately  the  same  amount  (expressed 
in  Nm3)  of  gas,  independently  of  it  is  filled  with  02  or  H2,  the 
resulting  benefit  is  the  possibility  to  supply  one  bottle  of  H2 
instead  of  4.5  bottles  of  02. 

Appendix  A.  Calculation  of  the  ideal  work  required  for 
oxygen  separation 

The  ideal  work  required  to  separate  oxygen  and  made  it  available 
at  pressure  pF in  in  an  isothermal  process  can  be  calculated  consid¬ 
ering  the  device  shown  in  Fig.  A.1(0):  a  piston  sliding  in  a  cylinder 
closed  by  an  ideal  membrane  (M)  perfectly  permeable  to  02  but 
impermeable  to  the  other  air  components.  Initially  the  piston  is 
at  the  top  dead  center  in  contact  with  the  membrane.  By  moving 
leftward  the  piston  (Fig.  A.l),  the  volume  in  the  cylinder  increases 
and,  since  the  membrane  is  in  contact  with  air  at  atmospheric  pres¬ 
sure,  the  volume  is  filled  by  the  02  flowing  through  membrane  at 
the  partial  pressure  of  02  in  air  (0.21  atm).  The  work  produced  by 
the  piston  in  a  reversible  process  is  l\  =  -  p02  ^1  (the  negative  sign 
means  that  work  is  applied  by  the  fluid  to  the  external  environment) 
where  V\  is  the  volume  swept  by  the  piston  in  a  complete  stroke. 
When  the  piston  reaches  the  bottom  dead  center,  the  cylinder  is 
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closed  with  a  lid  that  prevent  the  02  inside  the  cylinder  to  cross  the 
membrane.  Then  the  piston  is  moved  rightward  to  compress  the 
fluid  to  pressure  pF in  (Fig.  A.l(2)). 

By  noting  that  the  process  is  isothermal,  compression  work 
required  in  this  phase  can  be  computed  as: 
fV2 

L2  =  -  pdV  (A.l) 

JvA 

Provided  that  the  fluid  behavior  can  be  described  by  the  ideal  gas 
law,  pV  =  Nq2RT ,  the  differential  dV can  be  expressed  as: 


dV  =  - 


Nq2RT 
P 2 


dp 


(A.2) 


where  Nq2  is  the  number  of  02  moles  contained  in  the  cylinder.  By 
substituting  (A.2)  into  Eq.  (A.l )  it  follows: 

pV2  pPFIN  j 

L2  =  -  pdV  =  N02RT  -E=NoRTlnEM  (A.3) 

Jvt  Jp  02  P  P°2 


In  the  last  phase  of  the  process,  the  pulsion  work  required  to 
eject  the  02  from  the  cylinder  is  simply  calculated  as  I3  =Pfin^2- 
The  overall  work  required  is  therefore  L  =  Li  +I2  +L3  but  since  the 
whole  process  is  isothermal,  I3  =  -L\  and  consequently  L  =  L2. 

Finally  specific  consumption  required  in  an  ideal  isothermal  pro¬ 
cess  that  provides  pure  02  atpFiN  is  given  by  the  expression: 

£rev=  2-  =i?rin5M  (A.4) 

no2  P  o2 


and  coincides  with  the  value  obtained  in  Eq.  (6). 


Appendix  B.  Calculation  of  pressure  losses  and  ejector 
model 


As  reported  in  Section  2.1 ,  due  to  the  high  operating  temperature 
(about  750  °C),  the  recirculation  of  the  H2-H20  mixture  from  the 
SOEC  cathode  outlet  to  the  SOFC  anode  inlet,  is  given  by  an  ejector, 
instead  of  a  traditional  blower.  A  preliminary  design  of  this  appara¬ 
tus  can  be  made  after  calculation  of  the  pressure  losses  within  the 
SOFC-SOEC  loop  (streams  6-8)  which  determines  the  conditions  of 
the  ejector  driving  flow  -  in  terms  of  flow  fraction  and  pressure  - 
necessary  to  provide  the  required  A P. 

Calculation  of  pressure  losses  is  carried  out  as  a  function  of 
the  flow  conditions  within  the  SOFC  and  SOEC  channels  which 
are  assumed  to  have  a  rectangular  frontal  size  of  100  mm  x  1  mm, 
according  to  [11  ].  Since  the  total  mass  flow  of  the  H2-H20  mixture 
is  fixed  by  the  value  of  the  desired  02  yield,  the  pertinent  mole  flow 
value  for  each  channel  and  the  number  of  channels  are  calculated, 
given  the  density  of  the  mixture  (function  of  p  and  T,  according  to 
the  ideal  gas  law),  assuming  an  inlet  velocity  of  0.5  ms-1.  The  length 
of  each  channels  can  consequently  be  obtained  from  the  total  cell 
active  area  which  is  calculated  by  the  ratio  between  the  current 
(which  is  also  fixed  by  the  desired  02  yield),  and  the  current  density 
of  the  stack. 

For  the  case  considered  in  Table  1  with  a  current  density  of 
500 Am-2  the  active  area  is  24.1  m2(g02  s-1)-1,  resulting  in  183 
channels  per  g02  s_1  with  length  of  132  mm. 

Once  the  channel  length  has  been  determined,  pressure  losses 
are  finally  calculated  following  a  procedure  already  reported  in  [13] 
based  on  the  fact  that  for  laminar  flow  in  a  gas  channel,  the  product 
of  the  skin-friction  coefficient  and  the  Reynolds  number  depends 
only  on  the  channel  dimensions  [14].  For  the  case  considered  a  pres¬ 
sure  loss  along  the  SOFC-SOEC  piping  of  about  420  Pa  has  been 
obtained. 

Given  the  pressure  losses,  it  is  possible  to  carry  out  a  preliminary 
design  of  the  ejector  system.  To  this  purpose  an  ejector  model  has 
been  developed.  With  reference  to  Fig.  B.l,  the  driving  or  primary 


flow  at  point  E  (representative  of  the  conditions  in  terms  of  total 
variables  of  stream  10  in  Fig.  6)  is  first  compressed  to  an  assigned 
pressure  by  means  of  a  compressor  with  0.7  isentropic  efficiency. 
The  flow  is  then  expanded  and  accelerated  in  a  convergent  nozzle  to 
the  pressure  of  the  secondary  flow  at  section  B.  The  secondary  flow, 
whose  conditions  in  terms  of  total  T  and  p  are  those  of  stream  6  of 
Fig.  6,  is  expanded  through  a  convergent  from  section  A  to  throat  B 
before  being  mixed  with  the  driving  flow.  The  two  flows  are  then 
mixed  (from  section  B  to  section  C)  and  the  resulting  flow  finally 
enters  the  diffuser  section  (from  C  to  D)  where  kinetic  energy  is 
converted  into  pressure  energy.  A  diffuser  and  convergent  efficiency 
of  0.87  and  0.97,  respectively  are  assumed  in  the  calculation.  These 
efficiencies  are  defined  as  follows: 


Pconv  — 


A /lre 
A  Jlis 


*7diff  = 


Ahjs 
A /lre 


(B.l) 

(B.2) 


where  A/ire  is  the  variation  of  the  static  enthalpy  along  the  actual 
transformation,  while  A/iis  is  the  variation  of  static  enthalpy  of  an 
isentropic  process  reaching  the  same  value  of  static  pressure  of  the 
real  case.  The  two  terms  A/iis  and  A /ire,  are  evidenced  in  Fig.  B.2  in 
case  of  real  and  isentropic  expansion  (line  AB  and  ABr,  respectively) 


s 

Fig.  B.2.  Expansion  process  through  the  convergent  A  and  B  (1).  The  variation  of 
static  enthalpy  of  the  real  process  (A hre)  and  of  the  isentropic  process  (A his)  is 
evidenced. 
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Fig.  B.3.  Results  of  the  ejector  model:  pressure  increase  as  a  function  of  the  velocity 
of  the  secondary  flow  for  different  sets  of  mass  fraction  and  pressure  of  the  driving 
flow. 

through  the  convergent  A  and  B  of  Fig.  B.l .  Analogous  considerations 
can  be  made  in  case  of  the  compression  through  the  diffuser  from 
C  to  D  (Fig.  B.2) 

The  model  assumptions  are: 

-  The  flow  is  one-dimensional 

-  The  mixing  process  occurs  at  constant  pressure 

-  Wall  shear  stress  are  neglected  in  the  throat  (from  section  B  to 
section  C) 

-  Nozzle  is  adiabatic 

The  model  input  data  are  the  following: 

-  mass  fraction  of  the  driving  flow  with  respect  to  the  total  flow 

-  total  temperature  and  total  pressure  of  driving  flow  (section  E) 
and  secondary  flow  (section  A) 

-  pressure  ratio  and  efficiency  of  the  driving  flow  blower 

-  convergent  and  diffuser  efficiency  defined  according  to  Eqs.  (B.l ) 
and  (B.2) 

-  velocity  of  the  secondary  flow  at  section  B 

The  model  applies  mass,  momentum  and  energy  and  conser¬ 
vation  equations  within  the  mixing  section  (B  and  C)  allowing  to 
determine  for  any  input  conditions  the  resulting  pressure  increase 


across  the  ejector  (from  section  A  to  section  D  in  Fig.  B.l).  Results 
of  this  analysis  are  reported  in  Fig.  B.3  where  the  pressure  increase 
is  plotted  as  a  function  of  the  velocity  of  the  secondary  flow  at  sec¬ 
tion  B  (which  is  representative  of  the  extent  of  the  expansion  of  the 
secondary  flow  through  A-B)  for  different  sets  of  mass  fraction  and 
pressure  of  the  driving  flow. 

As  expected,  the  curves  present  a  maximum  given  that,  as  the 
velocity  of  the  secondary  flow  on  section  B  increases,  in  the  rising 
branch  of  the  curve  entropy  losses  related  to  mixing  of  streams  at 
different  velocity  (from  section  B  to  section  C)  decrease  faster  than 
losses  increase  due  to  expansion  (in  primary  and  secondary  flow 
nozzle)  and  diffusion  (from  C  to  D).  It  can  also  be  noted  that  an 
optimized  ejector  can  provide  a  420  Pa  pressure  increase  (the  one 
required  to  overcame  the  pressure  losses  of  the  SOFC-SOEC  loop) 
by  adopting  a  driving  flow  mass  fraction  of  5%  and  a  1.8  pressure 
ratio.  In  this  case  the  velocity  of  the  secondary  flow  in  the  throat  is 
about  230  m  s-1  and  the  resulting  geometry  shows  a  proportion  of 
the  diameter  of  sections  A:B:D  equal  to  1.0:0.15:1.02  with  a  overall 
length  about  seven  times  the  inlet  diameter. 
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